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Abstract 



O . By a choice of new variables the pressure isotropy condition for 

spherically symmetric static perfect fluid spacetimes can be made a 
^jrj' quadratic algebraic equation in one of the two functions appearing 

in it. Using the other variable as a generating function, the pressure 
and the density of the fluid can be expressed algebraically by the 
^ ' function and its derivatives. One of the functions in the metric can 

■ also be expressed similarly, but to obtain the other function, related to 

the redshift factor, one has to perform an integral. Conditions on the 
generating function ensuring regularity and physicality near the center 
are investiagted. Two everywhere physically well behaving example 
solutions are generated, one representing a compact fluid body with a 
zero pressure surface, the other an infinite sphere. 



1 Introduction 

The aim of this paper is to discuss an algorithm that can be used to generate 
any number of physically realistic pressure and density profiles for spherical 
perfect fluid distributions without calculating integrals. Our key step is the 
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transformation of the field equation into a form algebraic in one variable. 
A similar equation was already written down in the paper of Burlankov[|l|. 
However, it was immediately transformed to isotropic coordinates, thereby 
obtaining a generating formalism similar to that of Glass and GoldmanP|0, 
where even the calculation of pressure and density requires integrals. 

Because of their inability to support shear stresses static perfect fluid 
configurations are expected to be spherically symmetric. Theorems show- 
ing this have been proved under fairly general conditions [^1^. In the static 
spherically symmetric case the only relation one gets from the Einstein equa- 
tions is the pressure isotropy condition. This condition together with a fluid 
equation of state and a central density value can be used to determine a 
spherically symmetric regular matter distribution, which either connects to 
an exterior vacuum region through a zero pressure hypersurface or extends 
to the whole spacetime. The uniqueness of this spacetime has been proved 
by Rendall and Schmidt for monotonic equation of states. Regularity of 
the metric can be shown even without the monotonicity assumption^. 

Since the addition of even the simplest equation of state makes the task 
of finding an exact solution extremely difficult, the usual method to proceed 
is to look for some solution of the pressure isotropy condition under some 
mathematical assumptions and hope that the resulting equation of state will 
be simple and physical. The most common way to choose the radial coor- 
dinate r is by setting the surface of the isometry spheres to Airr'^. Then by 
appropriately choosing the form of the two functions describing the metric, 
the pressure isotropy condition is a first order linear equation in one of the 
two variables. As already was pointed out by Wymanp|, for arbitrary choice 
of the other function this equation can always be solved by quadratures. In 
the other variable the pressure isotropy condition is either a second order 
linear differential equation or a first order Riccati equation, depending on 
the exact choice of the function. Considering this relative mathematical sim- 
plicity, it is not surprising that more than hundred static perfect fluid exact 
solutions have already appeared in the literature [^]. Unfortunately a large 
part of these solutions are unphysical in several aspects. Many do not have 
a regular center, others have no positive pressure and density, violate the 
dominant energy condition, or have unphysical sound speed. 

Even though the fleld equation is in principle solvable by quadratures, the 
resulting integrals can be calculated in terms of elementary functions only in 
some exceptional cases. The imposition of physicality near the center makes 
the possible generating functions even more complicated and the chance of 
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being able to calculate the integrals even less. There are some techniques 
that can be used to generate new solutions from known ones published in 



the literature |]T0[ |Tl| , but they also require calculation of integrals or solution 



of differential equations. Using isotropic coordinates, Kuchowicz[|l2[ writes 
the pressure isotropy condition into a form which is algebraic in one of its 
two variables. Glass and Goldman]^ introduce pressure and density re- 
lated variables, also in isotropic coordinates, to obtain a similar algebraic 
equation. They use one of the two variables as a generating function, and 
obtain formulae for the pressure, the density and the metric functions in- 
volving integrals of functional expressions of the generating function and its 
derivatives. They also give conditions on the generating function which en- 
sure that the resulting metric is physically realistic near the origin. In the 
present paper we follow a similar approach in area coordinates r. The main 
advantage of our proposed method is that it does not require the calculation 
of any integrals for the expression of the pressure, the density and one of the 
metric functions. A similar equation to our main equation has been already 
calculated by Burlankov||I[], but a transformed form of it was used there to 
generate solutions. 

In Section |^ we transform the pressure isotropy condition into a quadratic 
algebraic equation in one of its variables and propose to use the other variable 
as a generating function. The pressure, the density and one of the functions 
in the metric are expressed as algebraic expressions of the generating func- 
tion and its first and second derivatives. The other function, determining the 
redshift factor, appears only in differentiated form in the field equation, and 
consequently only its derivative can be given in an algebraic form. The func- 
tional form of the generating function is calculated for some known solutions 
in Section 0. A power series expansion near the center is used in Section ^ in 
order to establish necessary conditions on the generating function to make the 
center regular, the density and pressure positive, having a local maximum, 
and satisfying the dominant energy condition near the center. The metric 
induced by the simplest physically realistic polynomial generating function 
is calculated in Section |. The solution represents a fluid sphere with a reg- 
ular center and a zero pressure surface. It has physical density and pressure 
and casual sound speed for some choice of the parameters. In Section ^ a 
second simple choice of the generating function is used to calculate a further 
solution, for which the matter extends to infinity. 
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2 Field equations 



The metric (yf^jy of a general stationary spherically symmetric configuration 
can be written in area coordinates as 

ds' = -e^'^de + \-dr^ + [d^^ + sin^ M^^) (1) 

where v and B are functions of the radial coordinate r. We assume that the 
spacetime region is filled with perfect fiuid, and the fiuid velocity vector is 
proportional to djdt. Since solutions of the Einstein equations with nonzero 
cosmological constant can be interpreted as solutions with negative pressure 
or density contributions, we do not include a cosmological term explicitly in 
our equations. Calculating the Einstein equations G^^y = StcT^^, we obtain 
the energy density, the radial and the angular directional pressures as 

(l-B- rB') , (2) 
{2rBv' + B-l) , (3) 
2rBv" + 2rB{v'f + rv'B' + 2Bv' + B'] , (4) 
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p^ = 


IGvrr 



where the prime denotes derivatives with respect to the radial coordinate r. 
The only field equation one gets is the pressure isotropy condition p = = 
o,9 which takes the form 



r 



rz/' + 1)B' + [2rV" + 2r^{iy'y - 2rz/' - 2] 5 + 2 = . (5) 



This equation together with an equation of state p = p{fi) and a central den- 
sity /ic determine a unique spherically symmetric perfect fiuid configuration, 
at least in the ^ > case0. In this paper, however, we try to find the 
general solution of the pressure isotropy condition pr = p^ and leave the task 
of interpreting the resulting equation of states as a second step. Equation 
d^) is first order and linear in B, and consequently, for any given function u, 
the general solution for B can be given in terms of integralsj^. We instead 
proceed by transforming (^) into a simpler form. 

Because of the freedom in constant rescaling of the time coordinate t, only 
the derivatives of the function u appears in the expressions for fi, pr and p^. 
This would suggest to introduce z/' as a new variable. However, considering 
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the coefficient of B' in it turns out to be more useful to introduce a new 
function as 

f3 = ru' + 1 . (6) 
Then the field equation takes the simple form 

r(3B' + 2rBl3' + 2/5^5 - 8/35 + 45 + 2 = , (7) 

and the radial pressure becomes 

Pr = ^ mB -B-1) . (8) 
Introducing a further new variable 

a = (3'^B , (9) 
equation (@) becomes a second order algebraic equation in /3, 

2(a + l)/5^ + (ra'-8a)/5 + 4a = . (10) 

If instead of a one introduces its square root, z = ^/a as a new function, 
and denote \/B = b, then it is possible to get from (|^) an equation algebraic 
in b, which is already given in the paper of Burlankov|I|], 

2b^ + {rz - 4:z)b + z"^ + 1 = . (11) 

Although at first sight this equation appears as a more comfortable form 
directly giving the metric function b = \^ from the generating function 
z = ^/a, it has serious disadvantages. First of all, for most known perfect 
fluid solutions B and a are not a square, and the appearance of the further 
square roots makes the calculations even more cumbersome. The more se- 
rious second problem is that the simplest polynomial choices for z do not 



appear to give appropriate results. For example, the quadratic equation (p!T[ ) 
has no real solution for b when choosing z = 1 — ar"^ with a constant, while 
a = 1 — ar"^ gives the Einstein static universe, as we will see in Section |^. 
In general, introducing any functional form of a, such as a^, ^ or e°, as a 
new variable would give a different algebraic equation in the other variable. 
This makes our formalism non-unique in one hand, but gives opportunities 
to generate even more variety of new solutions on the other hand. 

A further disadvantage of the equation ( p^l]) compared to (|T0|) that it is 
nonlinear in both variables. Since equation (p!0D is linear in a, if the function 
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V and consequently /3 = rv' + 1 is given, one can express its general solution 
for a using integrals, 



a 



e''f> (c -2 j ^e'^dr^ , (12) 

= J^{f3'-A(3 + 2)dr, (13) 

where C is some constant. Unfortunately, there is no guarantee, that for 
those u for which the integrals can be given in terms of elementary functions, 
the resulting equations of states will be also physically reasonable. Instead 
of trying to investigate this further we focus on the alternative approach, 
considering the function a as the basic quantity. 
For any function a for which 

(8a - ra')^ > 32a(a + 1) (14) 

the quadratic equation (^) has two solutions for the function /?. We denote 
them by /?+ and 

— ra ± \J (8a — ra')^ — 32a(a + 1) 



p± = — 8a - ra' ± J {8a - ra'Y - 32a(a + 1) . (15) 

4(a + 1) 

The condition (|T^ will turn out to be not particularly restrictive, since as 



we will see in Section ^, at a regular center a — > 1 and then (|1^) holds as 
an equality at r = 0. Furthermore, we will see that positivity of the fluid 
density will ensure the existence of the square root near the center. Actually, 
it will also turn out that the root /?_ always belong to non-positive densities, 
and hence it is unphysical. For actual calculations it may be convenient to 
use the identity 

P^P- = ^ (16) 
a + 1 

to eliminate square roots from the denominator. 

Using (^ and (^, the metric functions -B+ and belonging to and 
and z/_ belonging to can be calculated as 



_ a _ (a + 1)' 2 



1 

I -(/3±-l)cir + C± , (18) 
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where C+ and C_ are constants determining the scahng of the time coordi- 
nate t. The integral (|1^) generally cannot be given in terms of elementary 
functions. However, as can be seen from (Q) and the physically impor- 
tant pressure and density can be calculated without performing integrals. 
Denoting the pressure and density belonging to /5+ by p+ and and those 
belonging to /3_ by p_ and 

P± = ^ (2/3±i?± - 5± - 1) , (19) 

/i± = ^(l-B^-rB'A . (20) 

The form of the equation of state is extremely complicated in general, and 
possibly cannot even be put into an r independent /x = /x(p) or /(/i, p) = 
form using elementary functions. However, the important point is that, in 
principle, all static spherically symmetric perfect fluid solutions, with all 
possible equations of states, could be given by suitably choosing the function 
a. Indeed, for any known solution, first f3, then a can be calculated in terms 
of the metric functions z/ and B, using the equations (^ and (P). 



3 Some known solutions 

In this section we look at the form of some known exact solutions in our 
formalism. Since the field equation (|l^) has in general two roots, most of 
these metrics are paired with a perfect fluid counterpart solution. 

The simplest example is the Minkowski spacetime with a = 1. This is 
the only constant a spacetime with a regular center. Then the quadratic 
equation ([T0|) has only one solution, P = 1. 

Next we consider the vacuum Schwarzschild spacetime with mass param- 
eter m. Then e^^ = B = 1 — 2m/ r, and from (^) and (^) 

^ = ^ . (21) 

r — 2m 

(m — r)^ m 3 r ^/ , s 

« = -f ^ = -^ + I-^-T7r^ + ^(^ ) • (22) 

r[r — 2m) 2r 4 8m lom^ 

However, given this form of a, ( pT] ) is only one of the solutions of (p!0|). The 



other solution, corresponding to the negative sign in ([15D is 



(r — m)(r — 2m) 

(3_ = 2^ '-^ 1 . (23) 

m^ — Amr + 2r"' 
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This belongs to a perfect fluid spacetime with pressure and density 



7jv'{lm — 4r) 

^" " 327rr3(r - 2mf ' ' 

m^(2r — 3m)(10r — 17m) . . 

^" " 327rr2(r - 2m)4 ' ^ ^ 

The solution has a complicated equation of state and does not have a regular 
center at r = 0. 

The next simplest solution with a regular center is the Einstein static 
universe with constant z/, with (3 = 1 and a = i? = 1 — ar^, where a is some 
positive constant. The other root of (|T^) is 

1 — ar^ , , 
= ^ . 26 

2 — ar^ 



Using ([T7|) and(|T8D, the functions in the metric (|l]) turn out to be 

(2 - ar^f 
4(1 -ar2) 



= C(2 - ar^) , i?_ = i^^I^l! . (27) 



The pressure 

a(3ar2 - 4) , , 

= 321(1-4 

can be positive for negative a, but the density 

a^r^(3ar^ — 5) , , 

= ^J^TT^ ^ 29 

327r(l — ar"^)"' 

is always negative near the center and vanishes for r = 0. Although this 
solution is unphysical in the zero cosmological constant case, it has a regular 
center, and the fluid obeys a remarkably simple equation of state 

87ra/i_ = -(a + 87rp_)(3a + 647rp_) . (30) 

The interior Schwarzschild solution is described by the metric components 



e^^ = h^a-\jl--\ , B = l--, (31) 
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where a, b and R are constants. The pressure is positive if 1 < a < 3. This 
metric can be obtained in our formahsm by choosing 



a 



i?2 



+ 



a - \ I 



1 + 



3 — a 



and taking the root 



n 2 



- 3a + 4 
4i?6(a - 1)3 



(32) 



^6 ^ Q(^^B 



r^-R^ + aR'^Jl- 



i?2 



(33) 



of the quadratic equation (pi)!). The other solution /5_ gives another perfect 
fluid metric with a regular center and a complicated equation of state. Un- 
fortunately, the matter density fi- turns out to be negative near the center 
when the pressure is positive. We will see in the next section that this is a 
general property. If one solution of (|10D is physically well behaving, then the 
other belongs to negative matter densities. 



4 Expansion around a regular center 

In order for the metric ([1|) to possess a regular center the functions u and B 
must have regular limits at r = 0. To obtain the standard area per radius 
ratio for small spheres the limit of B must be 1. Further restrictions may come 
from differentiability conditions in a coordinate system which is regular at 
the center, and also from the regularity of the pressure and density. It can 
be seen from (H) and (^ that the limit of the functions P and a at r = 
must be also 1. 

We take the expansion of the function a in the form 

oo 

a = J2 (^(f' , (34) 

where are constant expansion coefficients, and ao = 1. If we assumed 
that V and 5, and consequently a and /5, are analytic at the center in a 
rectangular coordinate system, then because of the spherical symmetry, all 
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their odd expansion coefficients would have to vanish, i. e. a, would have 
to be zero for odd i. We proceed without this assumption now, and examine 
the first few expansion coefficients one by one. 

If ai is nonzero, calculating (3 from ([T5|), B from ([T7|), p and fi from 
([l9| ) and (pO|), the leading term in the pressures p+ and p- turns out to 
be proportional to and the leading term in the densities /i+ and yU_ is 

proportional to If ai = then the pressure is regular but the densities 

are proportional to This shows that ai and 0:3 must be zero for regular 
matter distributions. No such restriction follows for and higher odd index 
coefficients. 

When «! = as = then the expression under the square root in ([1^) has 
the expansion 



(8a - ra'Y - 32a(a + 1) 
-48a5r^ 



= 4(a^-8a4)r^ 
16(a2a4 + 4a6)r*^ + 0(r^) 



(35) 



In order for f3 to exist, this must be non-negative, which holds near the center 
only if 8a4 < There are two different spacetimes belonging to a given a. 
One belongs to the positive sign in ([ISD , and the other to the negative sign. 
The expansions of the pressures belonging to /5+ and /3_ are 



87rp+ 



a2 



+ 



as 



a2 + a2 
3a2 



«2 



8a4 — 12az 



— 8a4 



87rp_ 




a2 a2 — a2-\/a2 — 8a4 — 12a4 I r 



3a2 



- 7 + 0(/ 



a2 — 8a4 



(36) 



(37) 



These readily show that the pressure can be positive only if a2 > 0. The 
expansions of the fluid densities are 



87r/i^ 



as 



ao 



5a4 — a2 



12a5 



a2 — 8a4 



(38) 



4a4 + 



a^ + 32a6 



72ai 



a2 — 8a4 (a2 — 8a4j 



r2 + 0(r3) 
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Since the positivity of the pressure requires a2 > 0, the density belonging 
to the negative sign in (|T3p, is necessarily negative near the center. This 
shows that if one does not consider spacetimes with a cosmological constant, 
the (3^ root of (|I^) is always unphysical. However, the density belonging 
to is positive if 04 < and 0:2 > 0. In this case the 804 < a;| condition 



automatically holds, and hence the expression under the square root in ([15D 
is always positive. The dominant energy condition fi > p also holds at the 
center if < —^a^. 

In a physical situation one expects the pressure and density to have a 
maximum at the center. This certainly holds for p+ in a neighborhood of the 
center, since the coefficient of in (PB]) is positive when 02 > and a4 < 0. 
If as > then /i+ also have a local maximum at the center. In the more 
reahstic case, when = 0, the decreasing nature of /i+ gives a condition on 

32aQ > - (aj- 404) y'ai - 8^4 - "2 > (40) 

which always holds when ag is positive. 

Functions a with as > seem to be physically realistic, except that their 
equation of state for the fluid has some bad properties. If an equation of 
state exists in the form p = p{fi) then 

dp dp du , ^ . 

— = — — (41) 

dr dp, dr 

must hold. However, at the center ^ = but ^ is nonzero. This shows that 

' dr dr 

^ must be zero when p takes its central value pc, which means zero sound 
speed there. It would be important to know whether a realistic monotone 
increasing equation of state p = p{p) rules out all odd index coefficients in 
the expansion of a. 
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5 A solution representing a compact fluid 
sphere 

The simplest physically realistic polynomial choice for the generating function 
a appears to be 

a = 1 + acr^ - (l - 2a^) , (42) 
8 ^ ' 

where a and c are positive constants satisfying < |, or rather < ^ in 
order to comply with the dominant energy condition near the center. For 
the sake of simplifying the square roots appearing from (|T5[) we introduce a 
new radial variable x defined by 

2 _ 4(sinx + a) 

c(l-2a2) ■ ^ ' 

Then the variable x is restricted by x > Xc = arcsin (—a) . The generating 
function a takes the form 

1-2 sin^ X , ^ ^ , 

" = ^3^. (44) 

It still contains the parameter c through 

sinx = (l - 2a^) - a . (45) 

Using that now 



r-^ = (sinz + a)-^ (46) 

dr cos X ax 



the function under the square root in (|T5|) becomes 

(8a-ra'? -32a(a + l) = (Scosx^^^^^^^) . (47) 

V 1 — 2a^ J 

The root of ( p!OD belonging to positive densities takes the simple form 

^ sinx + cosx 

(3+ = . 48 

cos X — a 

Using equations ( p!7D and ([T8|), one of the functions in the metric is 

^ (cosx-a)^cos(2x) 
+ (1 - 2a2) [1 + sin(2x)] ' ^ 

12 



while the derivative of the other simphfies to 



du^ cos X 

dx 2 (cos X — a) 



(50) 



This can be integrated to yield 



X a 



arc tanh 



1 + a X 
, tan — 



(51) 



where an integration constant has been absorbed into the scaling of the time 
coordinate t. 

Using X as the radial coordinate the metric takes the form 



ds" 



— exp 



X + 



2a 



I arc tanh 



1 + a X 
= tan — 



dt' 



[1 + sin(2x)] cos^ X 



c cos(2x) (sin x + a) (cos x — a)^ 

4 fsin X + a) / 2 . ■ 2 a-, 2\ 
^ {d'&^ + sm^ ^(p^j 



dx" 



(52) 



c(l-2a2) 



This metric appears to be a new solution. Since the scale of the time co- 
ordinate t is arbitrary, the constant c corresponds to a constant conformal 
transformation of the metric. The pressure and density are 



87rp 



= c 



(a — 3 sin x) cos x + (3a — sin x) sin x 

4(sina; + cosx) 
7 sin X — 11a + 2asin^a; — lOsin^x 
4(sina; + cosa;)^ 
3 + 2a^ — 4a sin ,x — 3 sin^ x + Qa siv? x 



(53) 



+c- 



2 cos a; (sin x + cos x) 



(54) 



The fluid pressure and density are monotonically decreasing out from a reg- 
ular center to a p = surface. The dominant energy condition fi > p is 



satisfied everywhere for a < y ^ ~ 0.514, and the sound speed ^ is positive 
and less than one for a < 0.184. Unfortunately it is very difficult to express 
the equation of state in /j, — ^{p) form, or even in an f{ii,p) — form. The 
simplest expression the author could obtain is a complicated polynomial, 
eighth order in p, a and c, and fourth order in /i. 
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6 An infinite gaseous sphere solution 



A further simple assumption on the generating function a is that it is the 
ratio of two polynomials of the radial coordinate r. Considering the results 
in Section ^, the lowest degree form which might give a physically interesting 
solution is 

a = 1 + ^lIL^ = 1 + a2 _ hr^ + - 6^^) + O(riO) , (55) 
1 + hr^ ^ ' 

where a and h are positive constants. It is convenient to introduce a further 
constant c defined by 

c^ = -{h-a^) (56) 

and use it in place of the constant h. The assumption that c is real and 
non- negative restricts the original constants into the range h > a?. In or- 



der to simplify the square roots appearing from ( |15|) while expressing /?, we 
introduce a new radial variable x defined by 

2 2c — Ssinha; 

^ " a2 (2 + c2)(2sinhx-c) ' ^ ' 

Then the center is at = arcsinhy, spatial infinity is at sinh ^ , 

and the new variable is restricted by < Xoo < x < Xc- The b < case 
could be treated in a similar way introducing sinus functions instead of sinus 
hyperbolics, but would lead to solutions which fail to satisfy the dominant 
energy condition /i > p at infinity. The generating function a takes the form 

4c + (c^ — 4) sinh x , . 

« = — ; ; • (58) 

(2 + c2)sinha; ^ ^ 

The function under the square root in ([15|) becomes a square 

,,o , , /8c (2c — Ssinhx) coshxX^ , , 

(8.-™r-32«(a+l)=(^p^-j^-^j . (69) 

The root of ( p!OD belonging to positive densities takes the form 

( c coth f — 2 ) cosh f 

/9+ = ^ , ' { .J . (60) 

cosh f + c smh | 
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Using equation (P!7[), one of the functions in the metric is 

[4c + (c^ — 4) sinhx] (^ctanh f + 1^ 



(2 + c2) (ccothf -2j sinhx 

while, from (|18D , the derivative of the other function becomes 

dvj^ c cosh X 

dx 4 sinh | (^cosh f + c sinh |) (c — 2 sinh x) 

This can be integrated to obtain 

1 



(61) 



(62) 



-In sinh— , 

2 2 2 (3 + c2) 

X 



. , / 2 + c tanh f 

2v4 + arc tanh ' ^ 



1 + c ) In cosh — he sinh — I — In (2 sinh x — c 



2 2. 

Using a; as the radial coordinate the metric takes the form 



(63) 



exp (2z/_|_) dt^ 



2c — 3 sinh x 



+ 



a? (2 + c^) (2 sinh x — c) 
cosh^ X 



(d^"^ + sin^ ^9ci¥?2) 
dx^ , (64) 



4a2 (2 + c2) (2c - 3 sinh x) (2 sinh x - cy B+ 

where is determined by equation (^) and Bj^ by (|6T|) . This second metric 
also appears to be a new solution. The constant a corresponds to a constant 
conformal transformation of the metric. The pressure and density are 



Slip 
Stt/z 



; (2 sinh X — c) (2 — 2c^ + 2 cosh x + 5c sinh x) 
4(ccosh f - 2 sinh |) cosh^ | 
a^(2 sinhx — c) 



(65) 



2 sinh I 



cosh^ t: cosh X 



3 7c^ - 12 sinh(3x) 



32 (^ccosh 

+6c {Sc^ - 4) cosh(3x) + 2 (12 + eic^ - 8c^) sinh(2x) 

+4c (62 - 13c^) cosh(2x) + (l56 - 475c2 + 32c^) sinh x (66) 



+2c (4 - 53c2) cosh X + 12c (l3c^ - 22 



As indicated by numerical plots of the quantities, the fluid pressure and 
density are monotonically decreasing to zero, out from a regular center to 
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infinity. For any choice of tlie parameters a and c, tlie dominant energy 
condition /i > p is satisfied, and the sound speed ^ is positive and less than 
one. Ehminating the variable x, one can obtain a complicated polynomial 
equation of state in the form f{^iiP) = 0, eighth order in p, fourth order in 
H, and twelfth order in a and c. Near spatial infinity, in the small ii and p 
limit, the equation of state is approximately linear, 

p 6 + 4v/4+^ 



7 Conclusions 

An algorithm has been given, which can be used to generate physically real- 
istic density and pressure distributions from a generating function a without 
calculating integrals. Any function a of which the first few expansion coef- 
ficients satisfy the simple conditions stated in Section ^ generate spacetimes 
which are physically well behaving at least in a neighborhood of the cen- 
ter. The resulting pressure and density distributions can contain arbitrarily 
many parameters, for example by choosing a to be a high order polynomial. 
Unfortunately, the actual forms of /x and p can be quite comphcated because 
of the square roots appearing. This makes the task of putting the resulting 
equation of state in a closed form very difficult. Unfortunately, a prescribed 
equation of state would yield a very complicated differential equation on the 
generating function. However, by trying many different functional forms for 
a, hopefully one could find configurations with simple and physical equations 
of states. 

The metric function v only appears in a differentiated form in the field 
equation and in the density and pressure expressions. Correspondingly, only 
its derivative v' can be given as an algebraic expression of a and its first 
and second derivatives. Because of this, in some sense, one could consider 
a generated solution as an "exact solution" even if only the derivative of v 
can be given in a closed form. However, for some simpler choices of a, the 
integral determining v can be calculated. This has been the case for the two 
example solution solutions presented in this paper. 

It would be important to find out whether the generating formalism could 
be made simpler, or the obtainable exact solutions and equations of states 
could become more physical, for example, by choosing a functional expression 
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f{a) of the generating function instead of a, or by using a different radial 
coordinate in place of the area coordinate r. 
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